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Abstract 

The dynamics of a probe D7-brane in an asymptotically AdS-Vaidya background has been investigated in the presence 
of an external magnetic field. Holographically, this is dual to the dynamical meson melting in the Af = 2 super Yang- 
Milles theory. If the final temperature of the system is large enough, the probe D7-brane will dynamically cross 
the horizon (black hole embedding). By turning on the external magnetic field and raising it sufficiently, the final 
embedding of the corresponding D7-brane changes to Minkowski embedding. In the field theory side, this means that 
the mesons which melt due to the raise in the temperature, will form bound states again by applying an external 
magnetic field. We will also show that the evolution of the system to its final equilibrium state is postponed due to 
the presence of the magnetic field. 
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I. INTRODUCTION 

A new phase of QCD, called quark-gluon plasma (QGP), is produced at RHIC and LHC by colliding two heavy 
nuclei such as gold or lead, respectively, at relativistic speeds [1]. Studying the properties of the QGP has attracted 
a lot of attention. For instance, the process of plasma formation and equilibration has been considerably investigated 
over the last couple of years. In fact it is generally agreed that the plasma thermalizes very rapidly. By thermalization 
one usually means that the plasma reaches the stage after which the hydrodynamic equations describe the evolution 
of the system, . Experimental observations indicate that such a plasma is strongly coupled and numerical simulations 
suggest that at the early stages of the plasma formation a strong magnetic field is produced [2]. The effect of the 
presence of a strong magnetic field on the properties of the plasma has been studied extensively in the literature [3, 4]. 

Strongly coupled nature of the QGP makes perturbative methods in studying its properties not applicable. A 
reasonable candidate to overcome this problem is AdS/CFT correspondence or more generally gauge/gravity duality. 
This duality states that a d-dimensional strongly coupled Yang-Mills theory is dual to a (d-l- l)-dimensional classical 
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gravity theory on asymptotically AdS space-time [5]. The corresponding field theory lives on the boundary of the 
dual theory (bulk gravity). 

In the AdS/CFT dictionary the vacuum and an arbitrary thermal state of the field theory correspond to pure AdS 
and AdS black hole, respectively. Thus the thermalization process in the field theory can be equivalent to the black 
hole formation in the bulk as has been shown in [6] . For instance the Vaidya metric which is the result of the collapse 
of the matter in the bulk can describe the dynamics of the plasma from zero to a finite temperature. 

The black hole formation in the bulk resembles thermalization in the gluon sector of the field theory. One can 
also study such a process in the meson sector, called meson melting. Depending on the temperature of the plasma, 
some of the quarkonium mesons (heavy mesons) stay stable and live in the plasma. The injection of enough energy 
can cause these mesons to melt. Holographically, mesons are described by the fluctuations of the shape of the probe 
D-branes in the asymptotically AdS backgrounds [7]. The mesons are stable if the probe brane does not touch the 
bulk horizon and they melt when it crosses the horizon [8]. The dynamical meson melting by which we mean the 
process of the meson melting in a time-dependent background (Vaidya metric) has been investigated in [9]. In this 
paper we are interested in studying the effect of a non-zero external magnetic field on the dynamics of the meson 
melting. We will also discuss how the magnetic field affects the equilibration time when the system gets to its final 
equilibrium situation. 


II. PROBE D7-BRANE IN THE VAIDYA BACKGROUND 

In this section we will study the set-up of the dynamical meson melting in the presence of a constant external 
magnetic field. According to the AdS/CFT dictionary the mesons in the gauge theory correspond to the fluctuations 
of the probe brane in the asymptotically AdS background. In fact, in the AdS-black hole background, the shape of 
the brane can be categorized into Minkowski Embedding (ME) and Black Hole Embedding (BE) which correspond to 
the mesonic (stable mesons) and melted (unstable mesons) phases, respectively [8]. The process of dynamical meson 
melting is described by the change in the shape of the probe brane dynamically. In order to achieve this goal we study 
the shape of the probe brane in the AdS-Vaidya background. It describes the black hole formation in the bulk which 
corresponds to thermalization in the gluon sector of the dual gauge theory. The background is time-dependent and 
therefore the shape of the brane in such a background can deform dynamically. 

Let us start with the AdS-Vaidya metric that is 

ds^ = GMNdx^dx^ = ^ [—F(V, z)dV‘^ — 2dVdz + dx^] + d(j)^ -\- cos^ (j) dfl^ + sin^ (j) dijP', (H-l) 


where 


F(V,z) = 1 - M(V)z^ (H.2) 

and we have set the radius of AdS space-time to be one, R = \. The above metric is written in Eddington-Finkelstein 
coordinate where the radial direction is represented by z and V shows the null direction. The boundary, where the 
gauge theory lives, is at z = 0 and V is the time coordinate on the boundary. M(V), which is an arbitrary function, 
represents the mass of the black hole which changes as time passes by until it reaches a constant value. The function 
for M(V) that we will work with in this paper, is 

r 0 V < 0, 

M{V)=Mf\ i[l-cos(f^)] O^Vs^AV, (H.3) 

[ 1 V > AV, 

where AV is the time interval in which the mass of the black hole increases from zero to Mf 
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that the radius of the event horizon is Xh = Mf. 

In order to add the fundamental matter in the gauge theory side we have to add the probe 
dynamics of the degrees of freedom living on the brane is explained by the DBI action as 

S = -Tr j d^^^/-gab + {2TTa')Fab , (II.4) 

where gab is the induced metric and Fab is the gauge field strength on the brane. a, b are the brane coordinates and 
gab is defined as 


which is constant. Note 
branes to the bulk. The 


gab = daX^dbX^GMN- 


(H.5) 
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We are following the same notations as ones introduced in [9] and will not repeat the details. The D7-brane is 
embedded along the six directions of the bulk metric, x and ^ 3 . We choose the other two coordinates on the brane 
to be null coordinates u and v. For the rest of the bulk coordinates we choose the ansatz 

V = V{u,v), z = Z{u,v), 4> = ^{u,v), '0 = 0. (II-6) 

Since we are interested in studying the effect of the external magnetic field on the dynamical meson melting, we 
choose the magnetic field to be 


F — R 

^ XxX2 — 1 


(II.7) 


that is constant. 

Since we would like the u and v coordinates to stay null, we have to impose the following constraint equations on 
the equations of motion for the fields on the brane, obtained from the DBI action 

= 0 , ( 11 . 8 a) 

FVl + 2Z^^Vu-Z^{Z-^)% = 0, (II. 8 b) 

where 4 '(m,u) = which, in fact, gives the shape of the brane. With the above assumptions the DBI action 

becomes 

/ cos^ $ I 1 

dudv^^^J ^ {FVuV, + . (II.9) 


The difference between the above action and the one in [9] is that we have considered non-zero magnetic field which 
appears in the square root. Therefore, the equations of motion read 


ZB~^ 3 

Vuv = + -tan(Z^-)[(Z^),„y„ + 


^ Fz- 


FBo 


V,uV„, 


FZS^ o S F 

Z,uv = -+ - tan(Z4-)[(Z«'),,Z„ + - -^VuV, 


— - ( F^z — 


FB, 


{FVuV, + VuZ,. + V(„Z„) 


'i’.uv = 


3tan(Z4') -^FB^ 




2Z 


^ FB 2 3tan(Z^')\ 
^ z -—I- 


, X , X Sr-I-3^^'tan(Z^') , X 1 




'SBt 

j {FVuV, + V^Z^^ + 


2 Z 


-v,y„ 


where 


= 1 ± 


S, =3 


1 -k (27ra')2^4^2 ’ 

2 

1 -k {2TTa'YZ~B^' 


(Il.lOa) 

(Il.lOb) 

(II.lOc) 

(Il.lla) 

(Il.llb) 


By setting B = 0 we can recover the results of [9]. 

The above equations of motion are solvable, in spite of being non-linear and time-dependent, if one uses numerical 
methods. Here we apply finite difference method as has been used in [9]. The solution to these equations will tell 
us how the shape of the brane changes due to the dynamical background and the presence of the non-zero magnetic 
field. According to the AdS/CFT dictionary the near boundary expansion of the shape of the brane gives the mass 
of quark (to) and the quark condensate (c). More specifically 


^iV,Z)\z^o = m{V) + 


c{V) + 



Z^ + .... 


( 11 . 12 ) 


Note that in the presence of the magnetic field we rescale all the parameters with to, for instance B is rescaled as 
This means that to = 1 in all our results. By solving the equations of motion we will be able to see how the quark 
condensate behaves in time. The boundary and initial conditions are summarized in the following subsection. 
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A. Boundary conditions and initial data 

In order to solve the equations of motion we need to define the appropriate boundary and initial conditions. To 
begin we have to distinguish between two configurations mentioned before; Minkowski embedding and Black Hole 
embedding. Since the ME does not cross the horizon, we have to impose two boundary conditions; one at the AdS 
boundary Z\u=v = 0, where we have used the same definition as [9], and the other at the pole where 

the 3-sphere on which the D7-brane wraps, shrinks to zero. For the BE, the 3-sphere never shrinks to zero as the 
D7-brane crosses the horizon and therefore we do not need to impose any boundary condition at the pole. In order 
to derive the appropriate boundary and initial conditions in the presence of the magnetic held, we follow the same 
path as [9] where we refer the reader to it for more details. 

• Boundary conditions at the AdS boundary 

The physical observables are obtained on the AdS boundary, where u = v and therefore we have to include 
the appropriate boundary conditions there. The boundary condition on Z and $ are simply Z\u^y = 0 and 
= m. In order to obtain the rest of the boundary conditions we expand the helds near the boundary and 
by imposing the regularity condition on the equations of motion and consistency with the constraint equations 
(II.8) we get 

— Vo(r') = 2Z ; Zuv\u=v = 0 I (II.13) 

dv 

where Vo('(^) = V\u=v Note that the boundary conditions are not affected by the presence of the external 
magnetic field. 

• Boundary conditions at the Pole 

Similar to the previous set of boundary conditions the presence of the magnetic field does not alter the boundary 
conditions at the pole, ^\u=v+^ = f for vanishing magnetic field. Therefore we have 

Z,u = , Vu = V^ , (11.14) 

a.t u = V + ^. 

• Initial data 

The background at which the probe brane is embedded is time-dependent where before the energy injection, 
E < 0, the space-time is pure AdS. Therefore, our initial data includes the static solution to the equations of 
motion for pure AdS with non-zero magnetic field [3, 4]. 


III. NUMERICAL RESULTS 

In this section we study the response of the system to the time-dependent change in the temperature in the presence 
of an external magnetic field. The response of the system is described by the behaviour of quark condensate, c, in 
terms of the boundary time. When the magnetic field is zero, this behaviour can be classified into three categories: 
Minkowski embedding. Black Hole embedding and Overeager case. 

• Minkowski Embedding 

For this case one starts from an initial ME at zero temperature and by increasing the temperature in the bulk 
the hnal configuration is still ME. Since the initial ME state represents the mass of the quark in the dual field 
theory, if the ratio of m to the final temperature of the system is much larger than one, we can be sure that 
the embedding will always stay Minkowski. Therefore by raising the temperature, the phase transition does not 
happen and the mesons survive in the system. In this case the quark condensate shows an oscillatory behaviour 
in time about its equilibrium value. 

• Black Hole Embedding 

If Y becomes smaller than one, even though the initial state is ME, the D7-brane crosses the horizon as the 
temperature is increased and thus the final configuration becomes BE. In this situation the quark condensate 
falls from the initial zero value to its hnal equilibrium amount. 
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FIG. 1: This figure (right) shows that how an overeager configuration crosses the horizon and jumps out of it and becomes ME. 
W(p) = z~^ sinc^ reperesents the shape of the brane [9]. Accordingly, the graph on the left shows that c{V) starts oscillating 
at the final stages, as it is expected from ME. This figure is plotted for B = 1.87140, AV = 0.5 and ru = 1.25. 




FIG. 2: c{V) for ME (left) and BE (right). In both figures blue solid curves show c{V) — Co, with B=0.56446, in which 
Co = 0.0747132. co is the initial configuration quark condensate, before the black hole formation in the bulk. The red dashed 
curves show c{V) in the absence of magnetic field. In all cases th = 1.25 and AV = 1.0. 


• Overeager Case 

Provided that the time-scale of the change in the temperature is small, if the corresponding mass of the initial 
ME configuration is larger but close to the final temperature, the overeager case happens. This means that the 
D7-brane crosses the horizon and returns back again to another ME. In this case c{V) decreases from zero to a 
minimum value and then starts increasing to its final equilibrium amount. 

Before reporting our results for the non-zero magnetic field we would like to mention an interesting observation in 
the overeager case. As can be seen in the figure 1, c{V) starts oscillating at the final stages of its time-dependent 
evolution. This oscillation is around a final equilibrium value which in fact is the c for a ME in the presence of the 
corresponding black hole in static case. This also confirms that the overeager solution ends up in ME. Looking at 
the plot of the shape of the brane, figure 1 (right), we can see that at the initial times after the energy injection the 
shape of the brane is described by the BE while at sufficiently later times it becomes ME. 

We have studied the effect of the presence of a non-zero external magnetic field on the above embeddings as 
presented in figure 2. Figure 2 left(right) compares the behaviour of the ME(BE) between the two cases. Note that 
the blue curves for non-zero B are shifted in order to have the starting point of c{V) plots at c = 0. As observed in 
the figure, the equilibrium value of c decreases for B ^ 0 and the oscillations happen with a smaller phase velocity. 

We have computed the power spectrum for the ME oscillations in both cases and present the results in table I. The 
power spectrum is obtained from the discrete Fourier transform of c{V). In fact the peaks of the power spectrnm 
give us the mass spectrum of the stable mesons. Interestingly we observe that the meson masses increases by raising 
the magnetic field as seen in table I. Since for the binding energy of the stable mesons we have Eb ^ y/XM^es [10], 
where A is the t’Hooft coupling constant, the binding energy increases in the presence of the magnetic field. We 
should emphasize that the results obtained here come form solving the non-linear equations of motion rather than 
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TABLE I: Bound-States or Stable Mesons obtained from the ME Power Spectrum 


Stable Modes 

B = 0 

B = 0.56446 


2.82052 

3.04731 

0J2 

4.83518 

5.18042 

UJ3 

6.84984 

7.21553 




V 


V 




V 


V 


FIG. 3: This figure shows how the BE changes to overeager and consequently to ME after one raises the magnetic field. We 
have set = 1.25 and AV = 0.5. 


considering the linear perturbation around the equilibrium shape of the brane. Figure 2 (right) confirms the fact that 
c{V) decreases by applying an external magnetic field. We will elaborate more on this in figure 4. 

A remarkable phenomenon has been observed in the presence of the external magnetic field. Let us assume we start 
from a ME which results in BE after the injection of energy when the magnetic field is zero, figure 3 up-left. If we 
raise the magnetic field, it can still give BE as can be seen in figure 3 up-right. Interestingly, by further increasing 
the value of the magnetic field, with the same initial configuration and time-scale of the energy injection, the final 
configuration becomes overeager, figure 3 down-left. If the magnetic field becomes large enough, the configuration 
stays ME and never touches the horizon, as shown in figure 3 down-right. It is instructive to see what this means in 
the field theory side. At zero magnetic field, due to the energy injection the mesons are melted; indicating that quarks 
and anti-quarks live freely in the plasma. Raising the external magnetic field causes them to form bound-states again 
and the mesons become stable. This is consistent with our physical intuition from the non-dynamical picture. As has 
been discussed in [4], this observation can be explained by the fact that due to the non-zero magnetic field the probe 
brane becomes resistant to change in its shape. In the field theory side this means that the binding energy of the 
mesons become larger and therefore, more energy is needed to melt the mesons. 

When we start from a ME which, due to the energy injection, ends up in a BE, there exists a time at which c{V) 
relaxes to its equilibrium value. We call this equilibration time Veq that the system is at its equilibrium state. We 
define a time-dependent parameter 


e{V) 


o(E) Ceq 

c{V) 


(III.l) 
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FIG. 4: The plot on the left shows c{V) for BEs for different values of B. The graph on the right shows e{V) with respect to 
the boundary time, obtained from (III.l), for the corresponding plots in the left figure. The equilibration times are represented 
in table II. We have set Vh = 1.25 and AV = 1.0. 


where the equilibration time is defined as the time which satisfies e{Veq) < 0.003 and e{V) stays below this limit 
afterwards. Note that Ceq can be obtained from the static solution at the corresponding magnetic field and final 
temperature. This function has been plotted in figure 4 for different values of the magnetic field. We can conclude 
from figure 4 (right) that the presence of the magnetic field delays the evolution of the system to its equilibrium state. 
The equilibration times corresponding to the figure 4 are given in table II which confirm, quantitatively, the previous 
statement. 


TABLE II: Equilibration Times with respect to the External Magnetic Field 


B 

Veq 

0.0 

3.27834 

0.3137 

3.33708 

0.5645 

3.49866 

0.7129 

4.8057 

0.7946 

5.15314 
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